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^ ; Abstract 
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ly^ I Exclusive semileptonic decays of B mesons to orbitally excited D mesons are 

considered beyond the infinitely heavy quark limit in the framework of the 
\^ • relativistic quark model based on the quasipotential approach. This model 

^p . agrees with the structure of heavy quark mass corrections predicted by the 

heavy quark effective theory and allows the determination of corresponding 
_i , leading and subleading Isgur-Wise functions. It is found that both relativistic 

Q^' and 1/mQ contributions significantly influence the decay rates. Thus, rela- 

Q^' tivistic transformations of the meson wave functions (Wigner rotation of the 

D . light quark spin) already contribute at leading order of the heavy quark ex- 

pansion and result in a suppression of S — > D^ev and B — > D\eu decay rates. 
On the other hand, the vanishing of the decay matrix elements at zero recoil 
/\ ' of a final D** meson in the infinitely heavy quark mass limit makes the l/mq 

C^ ' corrections to be very important, and their account results in a substantial 

enhancement oi B ^ Diev and B -^ D^ev decay rates. 
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I. INTRODUCTION 

The investigation of semileptonic decays of B mesons to excited D meson states is an 
important problem for heavy flavour physics. In particular, these decays can provide an 
additional source of information for the determination of the Cabibbo-Kobayashi-Maskawa 
matrix element Vcb as well as on the relativistic quark dynamics inside heavy-light mesons. 
The experimental data on these decays are becoming available now ||l]-0, and the B fac- 
tories will provide more accurate and comprehensive data. It is necessary to note that the 
experiment shows that only approximately 60% of the inclusive semileptonic B decay rate 
is due to the decays to ground state pseudoscalar D and vector D* mesons. Thus the rest 
of these decays should go to excited D meson and continuum states. 

The presence of the heavy quark in the initial and final meson states in these decays 
considerably simplifies their theoretical description. A good starting point for this analysis 
is the infinitely heavy quark limit, rriQ — > cxd |^. In this limit the heavy quark symmetry 
arises, which strongly reduces the number of independent weak form factors 0. The heavy 
quark mass and spin then decouple and all meson properties are determined by light-quark 
degrees of freedom alone. As a result the heavy quark degeneracy of energy levels emerges. 
The spin Sq of the light quark couples with its orbital momentum I {j = I ± Sg), resulting 
for P-wave mesons in two degenerate j = 3/2 states {J^ = l+,2+)[] and two degenerate 
j = 1/2 states (0+, 1+). The heavy quark symmetry also predicts that the form factors for 
B — > D**ev decays, where D** is a generic P-wave D meson stateQ, can be expressed in terms 
of two independent Isgur-Wise functions Pj. However, in the infinitely heavy quark limit 
the decay matrix elements between a B meson and an orbitally excited D meson vanish at 
zero recoil because of the heavy quark spin- flavour symmetry 0. The kinematically allowed 
range for these decays is not broad. As a result the role of relativistic and finite heavy quark 
mass contributions not vanishing at zero recoil is considerably more important here than in 
the decays to ground state D mesons. Thus, the magnitude of such corrections might be 
comparable with the leading order result. 

Recently the first order l/mg corrections to the exclusive semileptonic B decays into 
excited charmed mesons were investigated within the heavy quark effective theory (HQET) 
0. The structure of the l/mq corrections to decay matrix elements, which follows from 
QCD and heavy quark symmetry, was determined. It was found that at the first order of 
heavy quark expansion the B ^ Di, B —y Dq and B -^ D^ matrix elements do not vanish at 
zero recoil and can be expressed at this kinematical point in terms of the leading Isgur-Wise 
functions. Away from the zero recoil point new subleading Isgur-Wise functions arise, which 
cannot be determined from symmetry considerations alone. Thus for their determination, 
additional model dependent assumptions are necessary. In Ref. ||^ an estimation of these 
functions, based on the non-relativistic quark model as well as on some additional assump- 
tions, was made. It is just here that we can apply the relativistic quark model to get a more 



^Here J = j i: 1/2 is the total angular momentum, and the superscript P denotes the meson 
parity. 

^For concrete P-wave meson states the standard notations Di, D2, D^ and Dq are used. 



consistent calculation of relevant Isgur-Wise functions and decay rates. 

Our relativistic quark model is based on the quasipotential approach in quantum field 
theory with a specific choice of the quark-antiquark interaction potential. It provides a 
consistent scheme for the calculation of all relativistic corrections at a given f ^/c^ order and 
allows for the heavy quark l/mq expansion. In preceding papers we applied this model to the 
calculation of the mass spectra of orbitally and radially excited states of heavy-light mesons 
1^, as well as to a description of weak decays of B mesons to ground state heavy and light 
mesons ^,0. The heavy quark expansion for the ground state heavy-to- heavy semileptonic 
transitions |TD] has been found to be in agreement with model-independent predictions of the 



HQET. We considered the exclusive semileptonic decays of B mesons to orbitally excited 
D mesons in the infinitely heavy quark limit in ^^ and found the important relativistic 
contribution to the leading Isgur-Wise functions arising from the relativistic transformation 
of the meson wave function. 

The paper is organized as follows. In Sec. II we briefly present the necessary HQET 
results on the B -^ D** transition matrix elements obtained in Ref. 0. In Sec. Ill we 
describe our relativistic quark model, putting special emphasis on the calculation of de- 
cay matrix elements and on the relativistic transformation of a meson wave function from 
the rest reference frame to the moving one. The heavy quark expansion for decay matrix 
elements is carried out up to the first order I/ttiq corrections and compared to model in- 
dependent HQET predictions in Sees. IV- VI. There we present our results for leading and 
subleading Isgur-Wise functions and compare predictions for decay rates with and without 
the l/mq corrections being taken into account. Taking account of l/mq corrections leads 
to a substantial enhancement oi B ^ Dieu and B -^ D^ev decay rates and gives better 
agreement between theoretical predictions and available experimental data. We also present 
the electron spectra for the considered decays and test the fulfilment of the Bjorken sum 
rule. Section VII contains our conclusions. 



II. HQET RESULTS FOR DECAY MATRIX ELEMENTS 

A. B ^ Dieu and B —>■ D^^eu decays 

The matrix elements of the vector and axial-vector currents between B mesons and Di 
or D^ mesons can be parametrized in the following way 

{D,{v\e)\-crh\B{v)) , 

, = Me ^ + Uv2V^ + fv^v ^){t ■ v), 

{D,{v',e)\-cri.b\B{v)) _^.^,.,,^., 
— - Uas e^f/3f , 

■^ymDiTUB 

y/mD*rnB 

where v {v') is the four-velocity of the B [D**) meson, e^ {e^'^) is a polarization vector 
(tensor) of the final vector (tensor) charmed meson, and the form factors fi and fcj are 



dimensionless functions oi w = v -v' . The double differential decay rates expressed in terms 
of the form factors read as follows |^J^ : 



d^r 



D, 



dwdcos 9 



d^r^. 

dwdcos 6 



SFor^v/^;;^^! sin^ 6 {w - n)fv, + {w' - l){fv, + nfv,) 



+ (1 - 2nw + rl) (1 + cos2^)[/^ + {w^ - l)/2] - 4cos^v/^;;^^T/yjA 



^'-Y^rliw'-lfl^ll^My^Q 



•4 
13 



[w - r2)kA^ + {w^ - 1)(/^A3 + r2k 



Ai) 



+ (1 - 2r2W + rl) (1 + cos^ e)[k\^ + {w^ - l)k^] -AcosOVw^ - IkA.ky , (2) 



where Fq = Cj^ \Vcb\^ m^ / {192tt^) , ri = rriD^/mB, T2 = 'm'£)*/'mB, and 9 is the angle between 
the charged lepton and the charmed meson in the rest frame of the virtual W boson. 

The main predictions of HQET for the structure of the B -^ Diev form factors look as 
follows [0: 

Ve/A = -{w + l)r - eb{{w - 1)[(A' + A)r - {2w + l)ri - t^] + (w + 1)17^} 

-£e[4(«;A' - A)r - ?>{w - l){n - r^) + (ti; + l)(r/ke - 2r/i - Sr/a)], 
Ve/y, = (1 - «;2)r - £5(«^' - 1)[(A' + A)r - {2w + l)ri - Ts + r]^] 

-£e[4(w + 1)M' - A)r - (w;2 - l)(3ri - 3r2 - 77ke + 2r/i + Sr/g)], 
Ve/y, = -3r - 3efe[(A' + A)r - {2w + l)ri - Ts + ^^ 

-£e[(4w; - l)ri + 5r2 + Sr^kc + lOr^i + 4(«; - 1)772 - 5r/3], 
Ve/vs = (^ - 2)r + €^{{2 + w;) [(A' + A)r - (2w + l)ri - r^] - (2 - w)r7t,} 

+ec[4(u;A' - A)r + (2 + w)ri + (2 + 3u;)r2 + (w - 2)r7ke 

-2(6 + w)^x - 4(w; - l)r/2 - (3«) - 2)773], (3) 

where Eq = l/{2mQ) and A(A') = M{M') — itlq is the difference between the heavy ground 
state (orbitally excited) meson and heavy quark masses in the limit ijiq -^ cx). The form fac- 
tor r is the leading order Isgur-Wise function (r is -\/3 times the function r3/2 of Refs. P,p!T[|). 
The subleading Isgur-Wise functions ti and T2 originate from the 1/mQ corrections to the 
b ^ c flavour changing current, while rike and rji form factors result from kinetic energy and 
chromomagnetic corrections to the HQET Lagrangian. 

The analogous formulae for B —>■ D2ei' have the form |^ 

kv = -T - EbliA' + A)r - {2w + l)ri - T2 + r]h] - edn - Ts + r/ke - 2r]i + 773), 
kA, = -(1 + w)t - eb{{w - 1)[(A' + A)r - {2w + l)ri - T2] + (1 + w)7]b} 



-Eciiw - l)(ri - Ta) + {w + l)(?7ke - '2r]i + 773)], 
kA, = T + EbiiA' + A)r - {2w + l)ri - ra + r/^] - edn + T2 - r^kc + 2r]i - 2r/2 - r^a). (4) 



B. B ^ DqCu and B ^ Dleu decays 

The matrix elements of the vector and axial currents between B mesons and Dq or D\ 
mesons can be parametrized as follows 
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where the form factors gi are functions of w. In terms of these form factors the double 
differential decay rates for B -^ D^eue and B — » D^eve decays can be expressed in the 
following way |^ 
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where Vq = Gp \Vcb\'^m% /{192tt^), Tq = m]:,*/mB and r^ = mo* / tub- 

The HQET predictions for the form factors of the decay B -^ D^ev are given by 
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5-- = C + ^c [Xke + 6x1 - 2{w + 1)X2] + £feXb- 
The analogous formulae for the decay B -^ D\ev look as follows 
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The form factor ( is the leading order Isgur-Wise function {( is twice the function T1/2 of 
Refs. [^|TI|). The subleading Isgur-Wise function (i originates from the l/rriQ corrections 
to the b ^ c flavour changing current, while Xke and Xi form factors result from kinetic 
energy and chromomagnetic corrections to the HQET Lagrangian. 

In the following sections we apply the relativistic quark model to the calculation of 
leading and subleading Isgur-Wise functions. 



III. RELATIVISTIC QUARK MODEL 

In the quasipotential approach, a meson is described by the wave function of the bound 



quark-ant iquark state, which satisfies the quasipotential equation |0 of the Schrodinger 
type m 
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where the relativistic reduced mass is 
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Here mq^q are the masses of light and heavy quarks, and p is their relative momentum. In 
the centre-of-mass system the relative momentum squared on mass shell reads 



b\M) 



[M^ - (mg + mQy][M^ - (m^ - ttiq)^] 
4M2 



(11) 



The kernel l^(p,q; M) in Eq. (^) is the quasipotential operator of the quark-antiquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, pro- 
jected onto the positive energy states. An important role in this construction is played by 
the Lorentz-structure of the confining quark-antiquark interaction in the meson. In con- 
structing the quasipotential of the quark-antiquark interaction we have assumed that the 
effective interaction is the sum of the usual one-gluon exchange term and the mixture of 
vector and scalar linear confining potentials. The quasipotential is then defined by [l^ 



V{p, q; M) = Ug{p)uQ{~p)V{p, q; M)uq{q)uQ{-q) 
= Uq{p)uQ{-p)l -as^^^(k)7^7Q 



+K!;nf(k)r^rQ;^ + v,u^) Mq)uQi-q), 



(12) 



where a^ is the QCD couphng constant, D^u is the gluon propagator in the Coulomb gauge 
and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 



u\p) 
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(13) 



>e(p)+m/ 



with e{p) = v^p2 + m'^. The effective long-range vector vertex is given by 



(14) 



where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the non-relativistic limit reduce 

to 
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^ ( 

conf \ 
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reproducing 



Konf(r) = Kif(r) + K 



' conf \ 



Ar + B, 



(16) 



where e is the mixing coefficient. 

The quasipotential for the heavy quarkonia, 



expanded in v'^/c^, can be found 



m 

Refs. |]I^ , pr5| and for heavy-light mesons in 0. All the parameters of our model, such as 
quark masses, parameters of the linear confining potential, mixing coefficient e and anoma- 
lous chromomagnetic quark moment k, were fixed from the analysis of heavy quarkonia 
masses fl^ and radiative decays [l^. The quark masses rrih = 4.88 GeV, rric = 1.55 GeV, 
rris = 0.50 GeV, rriu^d = 0.33 GeV and the parameters of the linear potential A = 0.18 GeV^ 
and B = —0.30 GeV have usual quark model values. The value of the vector-scalar mixing 
coefficient e = —1 has been determined by considering the heavy quark expansion |Ty] and 
meson radiative decays [16|. Finally, the universal Pauli interaction constant k, = —1 has 
been fixed from the analysis of the fine splitting of heavy quarkonia ^Pj- states |T^ . Note 
that the long-range magnetic contribution to the potential in our model is proportional to 
(1 + k) and thus vanishes for the chosen value of /t = — 1. 

In order to calculate the exclusive semileptonic decay rate of the B meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current J^ = C7^(l — 7^)6 
between a B meson and an orbitally excited D** meson takes the form |]17| 



(D"|j7(0)|i?)=/^^t^,.(p)r,(p,q)^B(q), 



(17) 



where r^(p, q) is the two-particle vertex function and "^3,0** are the meson wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
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frame. The contributions to F come from Figs. 1 and 2.^ In the heavy quark hmit rrib^c -^ oo 
only F*^^-* contributes, while F*^^^ contributes at l/mg order. They look like 

FW(p, q) = u,{p,h,{l - 7^)w,(g,)(27r)35(p, - q,), (18) 

and 

+^^P^ - q'^) g^(^) + gjg,) ^Q^9/-(^ - ^Q)}M(ib)uM, (19) 

where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, Q = c or 6, k = 
Pc - A; k' = q;, + A; A = p^.. - ps; e{p) = \w? + p2)V2; 



k^'\p) 



e{p) - (m7° + 7°(7p)) 



2e{p) 



Here [|T7| 



3 

Pc,q = ecMjf^ ± Y. n^'\vD")p\ 

Mb i=i 
and n*^*^ are three four-vectors given by 

The wave function of a P-wave D** meson at rest is given by 

^o- (p) = H\-) (p) = yf'M^) (p) , (20) 

where J and M are the total meson angular momentum and its projection, while j is the 
light quark angular momentum; ipD{j) (p) is the radial part of the wave function, which has 
been determined by the numerical solution of Eq. (^ in [0. The spin-angular momentum 
part 3^/*^ has the following form 



J 



xY,''-^^-''^XQ{^Q)x,{cr,). (21) 



^ The contribution F^^^ is the consequence of the projection onto the positive-energy states. Note 
that the form of the relativistic corrections resulting from the vertex function F*-^-* is exphcitly 
dependent on the Lorentz structure of the gg-interaction. 



Here {jimi, J2m2\J M) are Clebsch-Gordan coefficients, Y/^ are spherical harmonics, and 
x(o") (where a = ±1/2) are spin wave functions : 

x(l/2) = 

It is important to note that the wave functions entering the weak current matrix element 
(p!7|) are not in the rest frame in general. For example, in the B meson rest frame, the D** 
meson is moving with the recoil momentum A. The wave function of the moving D** meson 
\l/i3.. A is connected with the D** wave function in the rest frame "^d** q = ^D(i) by the 




transformation 17 



^D" a(p) = Dl/\RYjDl/\RYj^n-o{p), (22) 

where R^ is the Wigner rotation, L^ is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D^^'^(R) in spinor representation is given by 

(J i)^i?«A) = S'\p^,,)S{A)S{p), (23) 

where 



s{p)=\M^(^+ ""^ 



2m \ e{p) + mj 

is the usual Lorentz transformation matrix of the four-spinor. For electroweak B meson 
decays to 5- wave final mesons such a transformation contributes at ffist order of the l/rriQ 
expansion, while for the decays to excited final mesons it gives a contribution already to the 
leading term due to the orthogonality of the initial and final meson wave functions. 

IV. LEADING AND SUBLEADING ISGUR-WISE FUNCTIONS 

Now we can perform the heavy quark expansion for the matrix elements of B decays 
to orbitally excited D mesons in the framework of our model and determine leading and 
subleading Isgur-Wise functions. We substitute the vertex functions F*^^^ and F^^^ given by 
Eqs. (|TBp and (|TUp in the decay matrix element (^) and take into account the wave function 
properties (pOD-(p2D. The resulting structure of this matrix element is rather complicated, 
because it is necessary to integrate both over d^p and d^q. The S function in expression 
(p!8|) permits us to perform one of these integrations and thus this contribution can be easily 
calculated. The calculation of the vertex function F*^^) contribution is more difficult. Here, 
instead of a 5 function, we have a complicated structure, containing the Qq interaction 
potential in the meson. However, we can expand this contribution in inverse powers of 
heavy (6, c) quark masses and then use the quasipotential equation in order to perform one 



^Note that the quark model definition of '^d*(i/2) i^i (|20|), (p]) differs from the HQET one pLB] 
by an overall minus sign. 



of the integrations in the current matrix element. We carry out the heavy quark expansion 
up to first order in l/mg. It is easy to see that the vertex function T^'^^ contributes aheady 
at the subleading order of the l/mq expansion. Then we compare the arising decay matrix 
elements with the form factor decompositions (|l]) and (j^) and determine the corresponding 
form factors. We find that, for the chosen values of our model parameters (the mixing 
coefficient of vector and scalar confining potential e = —1 and the Pauli constant n = —1), 
the resulting structure at leading and subleading order in l/mq coincides with the model- 
independent predictions of HQET given by Eqs. (H), (^, (^, and (|). We get the following 
expressions for leading and subleading Isgur-Wise functions: 
i) B ^ Dicv and B ^ D2ev decays 



t{w) 



3 {W + 1)3/2 ] (2,,)3^^(3/2) \^ + Moimi^ + 1) 
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A' + A 

w + 1 
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w + 1 



{A' + A)t{w). 



(24) 

(25) 
(26) 



n) B -^ Dlei> and B —> Dleu decays 
y/2 1 



C(^) 



3 (u7 + 1)1/2 y (27r)3 
^ 2p 



d^p - ( 

^D(l/2) P + 
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MD(l/2)iw + l] 
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''dp 



eg + nig 



V^s(p), 



Ci(^) 



A* + A 



w 



aw) 



(27) 
(28) 



The contributions of all other subleading form factors, r]i{w) and Xii'^)^ to decay matrix 
elements are suppressed by an additional power of the ratio (w — l)/(w + 1), which is equal 
to zero ai w = 1 and less than 1/6 at tfmax = (1 + r^)/(2r) (r = ri, r2, Tq, or rl respectively). 
Since the main contribution to the decay rate comes from the values of form factors close 
to w = 1, these form factors turn out to be unimportant. This result is in agreement with 
the HQET-motivated considerations ||^ that the functions parametrizing the time-ordered 
products of the chromomagnetic term in the HQET Lagrangian with the leading order 
currents should be small. 

indicates that the derivative acts on the wave 



The arrow over d/dp in (p^ and 
function of the D** meson. All the wave functions and meson masses have been obtained in 
by the numerical solution of the quasipotential equation. We use the following values for 
HQET parameters A = 0.51 GeV, A' = 0.80 GeV, and A* = 0.89 GeV 0. 

The last terms in the square brackets of the expressions for the leading order Isgur- 
Wise functions t{w) ( P^ and C(^) (0) result from the wave function transformation ( ^2] ) 
associated with the relativistic rotation of the light quark spin (Wigner rotation) in passing 
to the moving reference frame. These terms are numerically important and lead to the 
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suppression of the ( form factor compared to r. Note that if we had apphed a simphfied non- 
relativistic quark model [§,|l^] these important contributions would be missing. Neglecting 
further the small difference between the wave functions ipD(i/2) and ipD{3/2), the following 
relation between r and ( would have been obtained 

C(^) = ^^(^)- (29) 

However, we see that this relation is violated if the relativistic transformation properties of 
the wave function are taken into account. At the point w = 1, where the initial B meson 
and final D** are at rest, we find instead the relation 

rW C(l) . 1 r ^V -^..,p)^^,^,p)_ (3„, 



V3 2 2 7 (27r)3^'' ' e, + m, 



obtained by assuming ipD{3/2) — V'd(i/2) — ^D"- The relation (^) coincides with the one 
found in Ref. |T^, where the Wigner rotation was also taken into account. 



V. NUMERICAL RESULTS AND PREDICTIONS 

In Table | we present our numerical results for the leading order Isgur-Wise functions r(l) 
and C(l) at zero recoil of the final D** meson, as well as their slopes pg^ = ~~^t and 

p^/2 = ~7^C 5 ill comparison with other model predictions P,p!9|-p4|. We see that most 
of the above approaches predict close values for the function r(l) and its slope P3/2, while 
the results for ({1) significantly differ from one another. This difference is a consequence of 
a specific treatment of the relativistic quark dynamics. Non-relativistic approaches predict 
C(l) ~ (2/v/3)r(l) (see (|2|)), while the relativistic treatment leads to {2/V^)t{1) > C(l) 
(see Eq. (pUD). The more relativistic the light quark in the heavy-light meson is, the more 
suppressed ( is with respect to r. We plot our results for leading (t(w), ({w)) and subleading 
{ti{w), T2{w), Cil"^)) Isgur-Wise functions for B -^ D**eiy in Figs. 3, 4 and for Bg — > D**eiy 
in Figs. 5, 6. 

We can now calculate the decay branching ratios by integrating double differential decay 
rates in Eqs. (0) and (^. Our results for decay rates both in the infinitely heavy quark limit 
and taking account of the first order l/mq corrections as well as their ratio 



R 



Bt{B -^ D**eu] 



mq- 



>oo 



are presented in Tables |1| and |ITI[ We see that the inclusion of l/mq corrections considerably 
influences the results and for some decays their contribution is as important as the leading 
order contribution. This is the consequence of the vanishing of the leading order contribution 
to the decay matrix elements due to the heavy quark spin-flavour symmetry at zero recoil 
of the final D** meson, while nothing prevents l/rnq corrections to contribute to the decay 
matrix element at this kinematical point. In fact, from Eqs. (|l|) and (^), we see that decay 
matrix elements at zero recoil are determined by form factors /vi(l), 5'+(l) and 5^^1(1), 
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Mii) = 


-8£e(A'-A)r(l) 


^+(1) = 


-^(^c + e.)(A*-A)C(l) 


9vA^) = 


(£,-3£fe)(A*-A)C(l). 



which receive non-vanishing contributions from first order heavy quark mass corrections. 
From Eqs. (|), ®, and (|) we find 

(31) 
(32) 
(33) 

Since the kinematically allowed range for these decays is not broad {1 < w < t^max ~ 1-32), 
the contribution to the decay rate of the rather small l/mq corrections is substantially 
increased. This is confirmed by numerical calculations. From Tables |T| and |TT| we see 
that the decay rate B — > D^^ev, for which all contributions vanish at zero recoil, is only 
slightly increased by subleading l/niq corrections. On the other hand, B -^ Dieu and 
B -^ DqCv decay rates receive large l/mg contributions. The situation is different for the 
B -^ Dleu decay. Here the l/rriQ contribution at zero recoil is not equal to zero, but it is 
suppressed by a very small factor (e^ — Se;,) (see Eq. (|33|) ), which is only ^ 0.03 GeV~^ for 



our model parameters. As a result the B -^ Dleu decay rate receives l/rnQ contributions 
comparable to those for the B -^ D^^ev rate. The above discussion shows that the sharp 
increase oi B ^ Dicv and B ^ Dq decay rates by first order l/mq corrections does not 
signal the breakdown of the heavy quark expansion, but is rather a result of the interplay of 
kinematical and dynamical effects. Thus we have good reasons to expect that higher order 
l/rriQ corrections will influence these decay rates at the level of 10 - 20 %. 

In Table H we present the experimental data from CLEO ^ and ALEPH p|, which are 
available only for the B -^ Dicv decay. For B -^ D2ei>, these experimental groups present 
only upper limits, which require the use of some additional assumptions about the hadronic 
branching ratios of the Dg meson. Our result for the branching ratio of the B — > Diei^ 
decay with the inclusion of l/ruQ corrections is in good agreement with both measurements. 
On the other hand, our branching ratio for the B ^> D2ev decay is only within the CLEO 
upper limit and disagrees with the ALEPH one. However, there are some reasons to expect 
that the ALEPH bound is too strong [||. 

In Table ^ we present our predictions for the ratios of decay rates B — > D2ev, B — »• 
DqCu, B -^ Dleu, and of the sum of all B — >■ D**ev decay rates to the rate B — > Dicv both 
in the limit mq —>■ cxd, and taking into account the l/rriQ corrections. It is reasonable to 



consider such ratios in order to normalize to a measured rate. In Ref. [|19| it is argued that 
a ratio Bt{B -^ D2eu)/BT{B -^ Dicv) = 1.55 ± 0.15 is a mere consequence of the heavy 
quark symmetry. In the heavy quark limit we confirm this result. However, the inclusion 
of l/mg corrections strongly influences this prediction and considerably reduce this ratio to 
a value close to 1. Such a reduction seems to be favoured by available experimental data. 
In the last row of Table |^ we give the sum of all B -^ D**ev branching ratios. We see 
that our model predicts that 1.45% of B meson decays go to the flrst orbitally excited D 
meson states. This result means that approximately 2.5% of B decays should go to higher 
excitations. 

In Fig. 7 we plot the electron spectra {l/Vo){dV /dy) for B -^ D**ez/ decays. Here y = 
2Ee/mB is the rescaled lepton energy. These differential decay rates can be easily obtained 
from double differential decay rates (H), (H), using the relation y = 1 — rw — r^w^ — 1 cos^ 
and then integrating in w over [(1 — y)^ + r^]/[2r(l — y)] < w < (1 + r^)/(2r). We present 
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our results both in the heavy quark hniit tuq -^ oo (dashed curves) and with the inclusion 
of first order l/niq corrections (solid curves). 



VI. BJORKEN SUM RULE 



Finally we test the fulfilment of the Bjorken sum rule |^3[ in our model. This sum rule 
states 

m m ^ 

where p^ is the slope of the B — > D^*^eu Isgur-Wise function, (^^^ and r*^™^ are the form 
factors describing the orbitally excited states discussed here and their radial excitations, and 
ellipses denote contributions from non-resonant channels. We see that the contribution of 
the lowest lying P-wave states implies the bound 

,= >UK(i)^^ 2^^.0.81, (35) 

4 4 3 ^ ^ 



which is in agreement with the slope p = 1.02 in our model |]I0[ and with experimental 



values p6|. 



VII. CONCLUSIONS 

In this paper we have applied the relativistic quark model to the consideration of semilep- 
tonic B decays to orbitally excited charmed mesons, in the leading and subleading order of 
the heavy quark expansion. We have found an interesting interplay of the relativistic and 
finite heavy quark mass contributions. In particular, it has been found that the Lorentz 
transformation properties of meson wave functions play an important role in the theoreti- 
cal description of these decays. Thus, the Wigner rotation of the light quark spin gives a 
significant contribution already at the leading order of the heavy quark expansion. This 
contribution considerably reduces the leading order Isgur-Wise function ( with respect to 
r. As a result, in this limit, the decay rates B — > DqCz/ and B -^ Dleu are approximately 
an order of magnitude smaller than the decay rates B -^ Dicv and B —>■ D2ev. On the 
other hand, inclusion of the first order I/itlq corrections also substantially influences the 
decay rates. This large effect of subleading heavy quark corrections is a consequence of 
the vanishing of the leading order contributions to the decay matrix elements due to heavy 
quark spin-flavour symmetry at the point of zero recoil of the final charmed meson. How- 
ever, the subleading order contributions to B ^ Dieu, B -^ DqCu and B -^ D^eu do 
not vanish at this kinematical point. Since the kinematical range for these decays is rather 
small, the role of these corrections is considerably increased. Their account results in an 
approximately twofold enhancement of the B — »• Diev and B —>■ D^ev decay rates, while 
the B —>■ D2ev and B — > D\ev rates are increased only slightly. The small influence of 
l/rriQ corrections on the B -^ D\ev decay rate is the consequence of the additional inter- 
play of l/rric and l/rrib corrections at the zero recoil point (see Eq. (p3D). We thus see that 
these subleading heavy quark corrections turn out to be very important and considerably 

12 



change the infinitely heavy quark hmit results. For example, the ratio of branching ratios 
Bt{B — > D2ei')/BT[B -^ Diev) changes from the value of about 1.6 in the heavy quark 
limit, rnq — > oo, to the value of about 1 after sub leading corrections are included. 

In conclusion, we have presented here the first self-consistent dynamical calculation of 
subleading heavy quark corrections in the framework of the relativistic quark model, which 
are found to be in agreement with the HQET predictions. 
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TABLES 

TABLE L The comparison of our results for the values of the leading Isgur-Wise functions r 
and (" at zero recoil of the final D** meson and their slopes py with other predictions. 





Ours 


i 


m 


m 


m 


0, n 


0], n 


r(l) 


0.85 


0.71 


0.97 


1.14 




1.02 


0.90 


pI/2 


1.53 


1.5 


2.3 


1.9 




1.5 


1.45 


C(i) 


0.59 


0.82 


0.18 


0.82 


0.70 ±0.16 


0.44 


0.12 


pI/2 


1.04 


1.0 


1.1 


1.4 


2.5 ±1.0 


0.83 


0.73 



TABLE IL Decay rates T (in units of |V;b/0.04|^ x 10"^^ GeV) and branching ratios BR (in %) 
for B -^ D**eu decays in the infinitely heavy quark mass limit and taking account of first order 
l/rriQ corrections, i? is a ratio of branching ratios taking account of l/w-g corrections to branching 
ratios in the infinitely heavy quark mass limit. 





mq 


— > oo 


With l/rnq 




Experiment 


Decay 


r 


Br 


V 


Br 


R 


Br (CLEO) gl Br (ALEPH) § 


B -^ Diev 


1.4 


0.32 


2.7 


0.63 


1.97 


0.56 ±0.13 ±0.08 ±0.04 0.74 ±0.16 


B -^ Dleu 


2.1 


0.51 


2.5 


0.59 


1.16 


< 0.8 < 0.2 


B -^ Dleu 


0.31 


0.073 


0.39 


0.09 


1.23 




B -^ D^eu 


0.25 


0.061 


0.59 


0.14 


2.3 





TABLE III. Decay rates T (in units of \Vcb/^M\'^ x 10"^^ GeV) and branching ratios BR (in 
%) for B -^ Dl*ev decays in the infinitely heavy quark mass limit and taking account of first order 
^/ruQ corrections. -R is a ratio of branching ratios taking account of 'i/ruQ corrections to branching 
ratios in the infinitely heavy quark mass limit. 







"iQ - 


-^ oo 






With l/niQ 






Decay 


r 






Br 


r 




Br 


R 


B -^ Dsiev 


1.5 






0.36 


4.5 




1.06 


2.9 


B ^ Dl^ev 


2.4 






0.56 


3.2 




0.75 


1.3 


B ^ Dl^eu 


0.53 






0.13 


0.77 




0.18 


1.4 


B ^ Dl^eu 


0.44 






0.10 


1.6 




0.37 


3.6 
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TABLE IV. Predictions for ratios of decay rates B -^ -D|ez^, B 



DqCu, and B 



Dfeu to the 



rate B -^ Diev in the mq -^ oo limit and taking account of l/mg corrections. In the last line we 
show our predictions for the sum of branching ratios (in %) of B decays to orbitally excited D** 
mesons. 



^Q 



CX) 



With l/mq 



Br(S -^ D*2ev)/Bx{B -^ Diev) 

Bt{B -^ D^eu)/Bi{B -^ Dieu) 

Bt{B -^ Dleu)/Bi{B -^ Dieu) 

EBr(S -^ D**eu)/Br{B -^ Dieu) 



1.59 
0.19 
0.23 
3.0 



0.94 
0.22 
0.14 
2.3 



J2Bi{B ^ D**eu) 



0.96 



1.45 
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FIG. 1. Lowest order vertex function F^^' contributing to the current matrix element (18). 
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FIG. 2. Vertex function F'^^ taking the quark interaction into account. Dashed lines correspond 
to the effective potential (p^. Bold lines denote the negative-energy part of the quark propagator. 



17 




FIG. 3. Isgur-Wise functions t{w) (upper curve), ti(w) (middle curve) and T2{w) (lower curve) 
for the B ^ Di2ev decay. 




FIG. 4. Isgur-Wise functions C(^) (upper curve) and Ci(^) (lower curve) for the B -^ D^^ev 
decay. 
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FIG. 5. Isgur-Wise functions t{w) (upper curve), ti(w) (middle curve) and T2{w) (lower curve) 
for the Bg — > Dgi 2^1^ decay. 
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FIG. 6. Isgur-Wise functions C(^) (upper curve) and Ci(^) (lower curve) for the Bs — > D*q ^ev 
decay. 
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FIG. 7. Electron spectra (1/ro) (dT/dy) for the B -^ D**eu decays. Dashed curves show the mq 
limit, solid curves include first order I/tuq corrections. 
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